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Abstract 

We propose a tomographic method to reconstruct the optical properties of a highly-scattering 
medium from incoherent acousto-optic measurements. The method is based on the solution to 
an inverse problem for the diffusion equation and makes use of the principle of interior control of 
boundary measurements by an external wave field. 
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The acousto-optic effect is a phenomenon in which the optical properties of a material 
medium are modified after interaction with acoustic radiation. Brillouin scattering from 
density fluctuations in a fluid l| or the ultrasonic modulation of multiply-scattered light in 
a random medium Q are familiar examples of this effect. It is well known that the scattered 
optical field carries information about the medium. This principle has been exploited to 
develop a hybrid imaging modality, known as acousto-optic imaging (AOI), which combines 
the spectroscopic sensitivity of optical methods with the spatial resolution of ultrasonic 
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171 ]. AOI holds great promise as 



a tool to probe hidden structure inside of highly-scattering media, such as clouds, paint 
and biological tissue — a problem which is of both fundamental interest and considerable 
applied importance [18 ] . For instance, in biomedical applications, optical methods provides 
unique capabilities to assess physiological function, such as blood volume and tissue oxy- 



genation 
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21] . At the same time, such methods have the molecular selectivity to map 



gene expression and other markers of bio-molecular activity. 

In a typical AOI experiment, a highly-scattering medium is illuminated by a coherent 
optical source and the resulting speckle pattern is registered by a detector. A focused ultra- 
sound beam is then introduced and the speckle modulation is recorded as the beam's focus 
is scanned throughout the medium. Since the scatterers in the medium are displaced by the 
acoustic wave, the scattered light undergoes a frequency shift which permits the localization 
of the resulting so-called tagged photons to the volume containing the focus. The intensity 
images that are obtained in this manner convey information about the medium. However, 
they are not tomographic, nor are they quantitatively related to the optical properties of 
the medium. 

In this Letter we consider the inverse scattering problem that arises in AOI. We show 
that it is possible to reconstruct tomographic images of the optical properties of a medium 
of interest from incoherent measurements of multiply-scattered light that is modulated by 
a standing acoustic wave. The principle advantages of the proposed method compared 
to conventional methods for imaging with diffuse light are two fold, (i) The resolution of 
reconstructed images is, in principle, much higher than in diffuse optical tomography (DOT). 
In particular, the inverse problem of AOI is well-posed and the image resolution is controlled 
by the acoustic wavelength. In contrast, the inverse scattering problem for diffuse waves is 
severely ill-posed, which leads to the relatively low resolution of DOT 221,123]. Physically, the 
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improvement in resolution in AOI can be understood to be a consequence of controlling an 
internal degree of freedom of the scattering medium (the density of scatterers) by means of 
an external wave field, (ii) Neither interferometric measurements of tagged photons nor the 
use of a focused acoustic wave field is required. This considerably simplifies the experimental 
realization of the method. 

We begin by developing a simple model for the acousto-optic effect in a random medium. 
The model accounts both for Brillouin scattering and multiple scattering of light due to 
density fluctuations and is formulated within the framework of radiative transport theory. 



That is, we do not account for coherent effects in multiple light scattering [241 ] . Consider the 
propagation of an acoustic wave in a fluid suspension of particles that scatter and absorb 
light. If the amplitude of the pressure wave is sufficiently small, each particle will oscillate 
about its local equilibrium position. We may thus regard the suspension as consisting of 
independent particles whose number density is spatially modulated due to the presence of 
the acoustic wave. If we focus our attention on a single neutrally buoyant spherical particle, 
its velocity u obeys the equation of motion 

"dt=V ?(v - u) - Vp ' (1) 

where we have not accounted for hydrodynamic interactions between the particles. Here p 
denotes the pressure, v is the velocity field in the fluid, a is the radius of the particle, p is 
its mass density, V = 47ra 3 /3 and r] is the viscosity of the suspension. Note that the neglect 
of hydrodynamic interactions means that, apart from a transient, the particle moves with 
the fluid. Consider a standing time-harmonic plane wave of frequency u with 

p = A cos(ut) cos(k • r + ip) , (2) 



where A is the amplitude of the wave, k is its wave vector and tp is the phase [25|. For 
simplicity, we have assumed that the speed of sound c s is constant with k = u/c s . The 
corresponding velocity field obeys the linearized Euler equation 

The oscillatory solution to JT]) is given by 

u = — sin(u;£) sin(k • r + ip)k , (4) 



and, as a consequence, we obtain for the position of the particle 

R = R 7; cosiut) sin(k • r + (p)k , (5) 

par 

where Ro is the equilibrium position. 

Let Ri, . . . , Rtv denote the positions of the particles in the fluid. Then, since each particle 
moves independently, we find that 

Rj = Roi 5- cos(a;t) sin(k • r + <z>)k , i = 1, . . . , N . (6) 

par 

The number density of particles is defined by 

g = Y,S(r-Mt)) ■ (7) 

i 

Expanding the delta-function to first order in the small parameter e = A cos(u)t) / (pc 2 ) and 
using ([2]), we find that g is given by 

g(r) = g (r) [1 + e cos(k • r + <p)\ , (8) 

where go(r) = J2i <K r — Roi) is the equilibrium number density of the particles. 

The propagation of multiply-scattered light is taken to be described by the diffusion 
approximation to the radiative transport equation. The electromagnetic energy density u 
obeys the diffusion equation 
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Dn 2 V 
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+ cp a u = in Q , (9) 



u + l— = 5(r-r ) on dfl , (10) 
on 

where n is the index of refraction of the fluid, p a is the absorption coefficient, I is the 
extrapolation length and r G dQ is the position of a unit-amplitude point source 0, 0]. 
The diffusion coefficient D is defined by 

D = — r t , (11) 

where p s is the scattering coefficient and g is the anisotropy of scattering. The scattering 
and absorption coefficients are related to the number density by p s = ga s and \i a = ga a , 
where a s and a a denote the scattering and absorption cross sections of the particles. We thus 
see that the optical properties of the medium are modulated by the presence of the acoustic 



wave. In addition, we account for variations in the index of refraction due to Brillouin 
scattering according to 

n(r) = n [1 + £7 cos(k • r + ip)] , (12) 
where uq is the index of refraction in the absence of the acoustic wave and 

7 = ~ (13) 

is the elasto-optical constant, with e being the dielectric constant of the fluid l||. Note that 
7 « 0.3 in water. Making use of the above definitions, we see that ([9]) becomes 

-V • D £ \/ip e + a^ e = in ft , (14) 

ifj e + £^± = g on 9ft , (15) 
on 

where ip e = u/n 2 , g = n 2 5(r — tq) and the modified absorption and diffusion coefficients are 
defined by a e = cn 2 fi a and D e = n 2 D. Using (jSJ) and (fT2"j) . we find that to first order in e 

a e (r) = a (r)[l + e(27 + l)cos(k-r + ^)] , (16) 
D e {r) = D (t) [1 + e(2 7 - 1) cos(k ■ r + cp)) , (17) 

where «o and D are the absorption and diffusion coefficients in the absence of the acoustic 
wave. The solution to (114"]) is given by 



^ e (r) = ^ (r) - e y dV[(2 7 + l)G(r, r')^o(r , )ao(r , ) 
+ (2 7 - 1) V r ,G(r, r') • V^o(r') A)(r')] cos(k • r + <p) + 0(e 2 ) , (18) 

where the Green's function G satisfies 

- V • £> (r) VG(r, r') + a (r)G(r, r') = 5(r - r') (19) 

and obeys homogeneous boundary conditions of the form (I1UI) with zero right hand side. We 
note that tf) e is proportional to the intensity that is measured by a point detector. 

We can estimate the magnitude of the acousto-optic signal from the above analysis. 
Consider a homogeneous volume of absorption a and diffusion coefficient D. It follows from 
(j2ip and fT22|) that the relative change in intensity due to the presence of an acoustic plane 
wave is of the order 

^ « e [l + (kL) 2 J e~ KL , (20) 



where L is the source-detector separation and k = ya/D. Choosing typical values of the 
above parameters in tissue: k — 1 cm -1 , L = 1 cm and e = 10 -3 , we find that AI/I « 10 -3 



which is expected to be observable [28]. We note AI/I would be significantly smaller for 
the case of a focused beam, assuming equivalent incident power. 

The inverse problem of AOI is to reconstruct ao and Dq from i/j e . To proceed, let us define 
= di/j/de\ e= o, which can be determined from measurements carried out in the presence 
and absence of the acoustic wave. Making use of (JIB"]) we obtain 

0(r) = J K(r, r') cos(k • r' + p)d 3 r' , (21) 

where 

K(r, r') = [(2 7 + l)G(r, r>o(i> (r') + (2 7 - l)V r ,G(r, r') • V^(r')i) (r')] . (22) 

Suppose we fix the positions of the optical source and detector and vary the wave vector k 
and the phase ip. Evidently, it is then possible to recover -^(r, r') by inversion of a Fourier 
transform. That is, 

d 3 k 

where the dependence of on k and (p has been made explicit. For simplicity, suppose that 
ao = 0. Noting that G and tpo depend upon Dq, it follows from (|22l that 

= ^[^o] , (24) 
where the nonlinear operator A is defined by 

A ^ - m - ijv^ ■ v^o • < 25 » 

Thus -Do is a fixed point of A, which can be found iteratively according to 

D^ = A[D^} , n = l,2,... , (26) 

where — > Dq as n — > oo, provided that ^4. is contracting. We note that at each step it 
is necessary to compute the Green's function G, which depends upon the current estimate 
of Dq. We further note that measurements from two independent sources are required to 
reconstruct both «o and D . 



K(r, r') = / -^e-^' [0(r; k, 0) + *0(r; k, 3vr/2)] , (23) 
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The above result provides an iterative solution to the inverse problem of AOL We now 
describe a direct method to solve the inverse problem. We begin by observing that <ft obeys 
the equation 



-V ■ A)V</> + [(2 7 - 1)V ■ A)W>o - (27 + l)cWo] cos(k ■ r + cp) + a o = in fl(J27) 

(f) + t d A = on dQ , (28) 
On 

which follows from ATI]) . (116]) and (TJI]). Next, we multiply ([27]) by ^0 and ([HD by 0, take 
the difference of the two equations that result and integrate over Q. After integrating by 
parts and employing the boundary conditions ( fTUl) and ( fT5l) . we obtain the identity 



(29) 



E(k,p)= / ^ (2 T -l) J D Q (VVo) 2 + (27 + 1)^0 cos(k-r + ^) 
where the surface term E is defined by the formula 

S(k, <p) = jD (r ) [0(r o ) + (2 7 - l)Vo(r ) cos(k ■ r + <p)] 

(27 - 1) 



/ d 2 rD ijj^ cos(k ■ r + cp) 
Jan 



(30) 



Since and ip are known on dQ from measurements, the function S can be determined. 
Provided that E is known for a sufficient number of values of k and <p, we observe that the 
Fourier transform of 



/(r) = (2 7 - l) J D (r)(V^o(r)) 2 + (2 7 + l)«o(r)^ 2 (r) 



(31) 



is known from available measurements. We first consider the case where ao = 0. We then 
find that 

A)(r) = M , , -o ■ (32) 



(2 7 -l)(V^ (r)) 2 
Making use of (jUJ), we see that ip obeys the nonlinear equation 



/ 



Wo) 



rV^o 



in 



/ j_ e d ^° an 
ip + I— — = g on oil . 

on 



(33) 
(34) 



This equation admits a solution if / is smooth and bounded from below by a positive 
constant. Once ipo is found by solving fl33|) . we can recover the diffusion coefficient D from 

©2D. 



Next we consider the general problem of recovering both «o and Dq. In this case, we 
require data from two sources g\ and g 2 , as specified by the boundary condition ffl5l) . We 
denote the corresponding solutions to (1T4"|) by ipi and ife. We also define 

/*(r) = (2 7 - l)D (r)(V^(r)) 2 + (2 7 + l)a (r)^ 2 (r) , k = 1,2 . (35) 

Solving for a and D we find 

, , = /i(r)(V^(r)) 2 -/ 2 (r)(V^(r)) 2 
01 J ~ (2 7 + l)[^(r)(V^(r))2-^ 2 (r)(V^i(r)) 2 ] ' 1 J 

01 J ~ (2 7 -l)[^(r)(V^(r)) 2 -^ 2 (r)(V^i(r)) 2 ] ' 1 J 

The ipk are then obtained by solving the system of nonlinear equations 

-V • D Vip k + a ipk = in tt , (38) 
4>k + i-g^ = gk on dQ, , (39) 

where a and D are defined by (|36|) and (|37|) . respectively and k = 1,2. Once the ^ are 
found, we can then recover a and D from fl36|) and fl37|) . which yields the solution to the 
nonlinear inverse problem of AOL 

Several remarks on the above results are necessary, (i) We observe that the inverse 
problem of AOI is well-posed. This can be seen from ff23l . which yields a bandlimited ap- 
proximation to the diffusion coefficient by an inverse Fourier transform. In contrast, the 
linear inverse problem of DOT is severely ill-posed, req uiring the inversion of a Laplace 
transform, a problem which has logarithmic stability ^l. I30I] . This ill-posedness is respon- 
sible for the relatively low resolution of images in DOT. Characterizing the stability of the 
inverse problem for AOI, along with developing numerical methods for (1381) . will be the 
subject of future research, (ii) The nonlinear equation (|33l) is of the form 

V-(/|V^r 2 V^)=0, (40) 



with p = 0. The case p = was also addressed in 3JJ in a different physical context and 



solved numerically using an iterative algorithm similar to the fixed point algorithm presented 



in (1261) . The case p = 1 was studied in 32L w here uniqueness results and numerical methods 



were described. See also related work in 



331 ] . (iii) The diffusion equation (Q is valid when 



the energy density varies slowly on the scale of the transport mean free path. This condition 



breaks down when the acoustic wavelength is sufficiently small. It would thus be of some 
interest to extend the theory we have developed to the transport regime. 

In conclusion, we have developed a tomographic method for acousto-optic imaging. Nei- 
ther interferometric measurements of tagged photons nor the use of a focused ultrasound 
beam is required. Our approach is based on the solution to an inverse problem for the 
diffusion equation with interior control of boundary measurements. 
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